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Abstract 

We consider a multicast scenario involving an ad hoc network of co-channel MIMO nodes in which 
a source node attempts to share a streaming message with all nodes in the network via some pre-defined 
multi-hop routing tree. The message is assumed to be broken down into packets, and the transmission 
is conducted over multiple frames. Each frame is divided into time slots, and each link in the routing 
tree is assigned one time slot in which to transmit its current packet. We present an algorithm for 
determining the number of time slots and the scheduling of the links in these time slots in order to 
optimize the connectivity of the network, which we define to be the probability that all links can achieve 
the required throughput. In addition to time multiplexing, the MIMO nodes also employ beamforming 
to manage interference when links are simultaneously active, and the beamformers are designed with 
the maximum connectivity metric in mind. The effects of outdated channel state information (CSI) are 
taken into account in both the scheduling and the beamforming designs. We also derive bounds on the 
network connectivity and sum transmit power in order to illustrate the impact of interference on network 
performance. Our simulation results demonstrate that the choice of the number of time slots is critical in 
optimizing network performance, and illustrate the significant advantage provided by multiple antennas 
in improving network connectivity. 
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I. Introduction 

A. Motivation 

Interference management in co-channel ad hoc networks is a challenging problem. Simple time-division 
multiple access (TDMA)-based designs are inefficient and usually result in relatively poor performance. 
While multiple antennas are useful in enhancing channel gain and reducing interference, incorporating 
the extra degrees of freedom offered by MIMO (multi-input, multi-output) nodes into the design of the 
network further complicates matters. Designs based purely on spatial-division multiple access (SDMA) 
are not appropriate for large networks, since the number of available antennas is usually insufficient for 
cancelling all of the co-channel interference. Consequently, space-time (STDMA) solutions must generally 
be employed, in which multiple network links are scheduled into each time slot, and beamforming 
techniques are used within each slot to mitigate the resulting interference. 

Prior work that has addressed STDMA scheduling for ad hoc networks has typically focused on finding 
solutions that maximize the sum throughput of the network iPTl-EI. However, such solutions inevitably 
lead to links with poor performance and localized network congestion, which cannot be tolerated in 
applications where the network must perform multicast streaming (5]. An alternative is of course to use 
techniques that ensure fairness (i.e., max-min rate, proportional-fair scheduling, etc.) O-Hl, but such 
techniques typically do not directly address network reliability. Performance may be fair, but how likely 
are the links able to achieve this performance? 

The goal of this paper is to suggest methods for addressing the above issues using physical (PHY) 
layer techniques in combination with interference-aware scheduling. We introduce a novel definition of 
network connectivity that quantifies the probability that all links in the ad hoc network are able to achieve a 
certain pre-specified throughput. The PHY-layer parameters (beamforming vectors, transmit powers) and 
the scheduling decisions are then made to maximize this connectivity metric, taking into account the 
interference produced by links that are simultaneously active. The scenario we have in mind is ad hoc 
multicast network streaming, in which a source node attempts to transmit a continuous data stream to all 
other nodes in the ad hoc network with maximum reliability. Applications of this problem include tactical 
military networks (source is the unit commander, other nodes are teams or individuals under its direction), 
sensor surveillance networks (source is a sensor streaming data from a detected event, network nodes are 
monitoring stations attempting to form a coordinated response), vehicular ad hoc networks (VANETs), 
etc. For example, a typical VANET scenario involves a vehicle in the network that detects an unsafe road 
condition that must be reported to all other network nodes. In such a case, throughput is not the most 
important aspect of the network, since the actual message may be relatively short. Instead, the reliability 
of the network in sharing the message from the source with all the vehicles in the network is the key 
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to achieving safety J9), ifTOl . Mobile sensing with distributed platforms (e.g., ground-based robots or 
UAVs) is another VANET application where connectivity is more critical than throughput, at least until 
a detection occurs. In sensing mode, such a network must continually be connected in order to arrive at 
consensus regarding the parameter(s) being sensed, or until one of them detects an object of interest, in 
which case the network may reconfigure itself for high-throughput data gathering. 

Connectivity is fundamentally different from both throughput and energy consumption, which are the 
metrics most commonly used to quantify the performance of a wireless network. Connectivity performance 
is more relevant for applications where robustness or reliability is the most critical factor, applications 
(e.g., such as in military or emergency response scenarios) where the overriding concern is ensuring that 
information is shared with all nodes in the network at some pre-determined minimum rate. Throughput 
and energy consumption do not address network robustness. Instead, we argue that the connectivity of 
the network, as defined in the paper, is a reasonable way to measure the network's robustness, and thus 
we set about choosing the network parameters to maximize the connectivity metric. A solution that 
achieved a higher throughput or lower energy consumption than our solution would necessarily have 
poorer connectivity and hence less robustness, and thus be a less desirable operating condition for the 
type of applications we assume. 

Our focus is on how the use of multiple antennas in the network can lead to dramatic improvements 
in connectivity. The additional spatial degrees of freedom offered by the antennas reduce interference 
between the links, they allow a "denser" scheduling of users for a given set of resources, and they reduce 
the amount of transmit power (and hence interference to other links) to achieve a given throughput. 
The combined gain of these effects can only be observed by considering the joint optimization of the 
transmit beamformers, scheduling and power control. Unlike individuals with hand-held communication 
devices, VANET nodes are typically not power constrained and are large enough to support the use of 
multiple antennas. Thus, vehicular communication systems are a natural platform on which to consider 
the performance of multiple-antenna communications. Since mobility is a defining aspect of VANETs, 
it is important to account for the impact of the resulting time-varying wireless channels. Unlike most 
work in multiple-antenna ad hoc networks, we explicitly take the time-varying channels into account and 
quantify their impact on the reliability (connectivity) of the network. 

B. Background 

This work draws on ideas and techniques that have been studied by many others, but in different con- 
texts, including connectivity, beamforming for interference networks, and interference-based scheduling. 
Relevant prior work in these areas is briefly discussed below. 



January 4, 2012 



DRAFT 



4 



Investigations of radio network connectivity have been conducted by researchers over the last several 
decades. The original connectivity paradigm was expressed in terms of the so-called "geometric disk" 
model and percolation theory. In the geometric disk model, two nodes are assumed to be directly connected 
if their distance is smaller than some minimum transmission radius. This results in a simple binary 
description of connectivity (i.e., the network either is or is not connected) but lacks an indication of the 
quality of information flow. Percolation theory revolves around finding node density conditions under 
which a given node belongs to an unbounded cluster of connected nodes iTTTTl . |[T2 1. However, neither of 
these approaches is reasonable for realistic networks, where one must consider the effects of fading and 
interference. The impact of fading on network connectivity has been addressed in |[lJl - lfT6*l . Of particular 
interest to this paper is lfl3l . which showed that multiple antennas can significantly reduce the node 
isolation probability in an interference-free network, and IPT61 , whose simulation results demonstrated 
that multiple antennas can significantly enhance the connectivity of an ad hoc network, measured as 
the number of links that meet a given requirement on the outage capacity, or the symbol error rate of 
orthogonal space-time block coding. Interference aspects of the problem have been studied in ifTTl , which 
investigated the connectivity of sensor networks with regular topologies. The network connectivity was 
defined as the probability that a path exits between any two pairs of nodes in the network, and simulation 
results illustrated how an increase in node density led to decreased connectivity due to interference effects. 

To reduce the impact of interference in ad hoc networks, interference graph and coloring techniques 
have been used in the design of scheduling or routing algorithms lTT8l - |[20ll . In lfl8ll . a linear programming 
(LP) method was proposed for computing lower and upper bounds for the maximum throughput that can 
be supported by a multi-hop network. A conflict (interference) graph was used to find the constraint 
conditions for the LP formulation. In [19], the authors proposed the construction of a link-directional 
interference graph to account for the directional traffic over each network link. They investigated a 
coloring algorithm with two colors on the interference graph to schedule transmissions in ad hoc networks 
employing TDMA or frequency division multiple access (FDMA). In [20} , active links in a multi-hop 
network are scheduled in an STDMA scheme where a frame is divided into equal length time slots and 
each time slot can be allocated to several links simultaneously. Utilizing the interference graph, a heuristic 
algorithm is proposed to minimize the frame length under a constraint that each link's minimum data 
rate requirement is satisfied. An earlier version of this paper [21] discusses the use of interference graphs 
together with MIMO for improved connectivity. 

The use of multiple antennas to improve the performance of ad hoc networks has been a topic of 
considerable recent research |[22l . For example, in |2"5lL the transmitter for each link uses the principal 
singular vector of the channel matrix as the transmit beamformer and each receiver uses MMSE beam- 
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forming to mitigate the inter-link interference, and it was shown that there exists an optimal number of 
active links that maximize the network throughput. Beyond this value, the network becomes interference- 
limited, and performance is degraded. The results in [24] show if each link transmits a single data stream 
without CSI, while the receiver uses partial-zero-forcing interference cancellation, the capacity lower 
bound increases linearly with the number of antennas. In (3J], an optimal scheduling policy is proposed 
that can maximize the average sum rate of the MIMO ad hoc network under the constraint that the 
average data rate of each link is larger than a certain threshold. In [7], two optimization problems are 
considered: one is to maximize the sum rate under sum power and proportional fairness constraints, and 
the other is to minimize the sum transmit power under a constraint on the minimum data rate of each 
link. In 12, several distributed scheduling methods were proposed for MIMO ad hoc networks. In these 
methods, the links whose channel condition satisfy a pre-defined threshold are divided into groups for 
simultaneous data transmission to maximize the overall network throughput. Also relevant is the recent 
research on MIMO interference networks, where techniques based on interference alignment [25] and 
game theory [26] have been proposed. 

C. General Approach and Contributions 

As mentioned earlier, in this paper we consider PHY-layer optimization and scheduling for an ad hoc 
MIMO network in multicast streaming mode. In particular, we assume a source node desires to stream 
data to all other nodes in the network via a pre-determined multi-hop routing table (for example, a 
minimum spanning tree). Our emphasis will be on obtaining high reliability and low congestion for the 
network by maximizing the network connectivity, which we define as the probability that all links are 
able to support the desired throughput. The contributions of this work are as follows: 

1) We examine the multi-hop multicast streaming problem using a detailed PHY-layer model, and 
demonstrate the significant benefit that multiple antennas, power control and proper scheduling can 
have on network robustness. 

2) We propose a new definition of ad hoc network connectivity that approximates the probability that 
all links in the network can achieve a certain average throughput. The metric provides a continuous 
measure of connectivity performance that is more descriptive than the simple binary metrics often 
used in wireless networks. In addition, this metric leads to a more robust solution for fast fading 
channels than techniques based on requiring that only the average rate of each link be above some 
threshold 0,0. 

3) Since optimal connectivity will require that network links be simultaneously active, we develop 
beamforming and power control algorithms for the MIMO interference channels in each time slot 
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that maximize the network connectivity metric. 

4) Together with the beamforming algorithm, we propose an approach for STDMA scheduling based 
on greedy coloring of the interference graph that finds near-optimal assignments of links to each 
time slot. 

5) Unlike nearly all prior work in this area, both the beamforming and scheduling algorithms take 
into account the fact that the channel state information (CSI) may be outdated when it is used, 
and we illustrate the impact of the outdated CSI on connectivity. This is particularly important for 
VANET applications, where the network nodes are mobile. 

6) We illustrate via a number of simulations the dramatic improvement in connectivity that can be 
obtained when the network nodes are equipped with multiple antennas. By jointly considering the 
problems of beamformer design, scheduling and power control, we observe that the use of multiple 
antennas provides a "multiplicative" benefit" that exceeds what one would expect from their use in 
addressing the problems individually. Furthermore, our simulations indicate that adding antennas 
to the network nodes actually reduces the relative performance loss due to outdated CSI. 

7) We derive analytic expressions for an upper bound on the network connectivity and a lower bound 
on the sum transmit power of the network assuming interference-free transmissions and no CSI 
errors. These optimistic bounds provide benchmarks to determine the robustness of the proposed 
approach. 

We note here that a key limitation of the proposed approach is that the optimization described above 
takes place assuming that the multicast routing tree is pre-determined and remains fixed. This is clearly 
suboptimal since the routing decisions determine the links, which in turn create the interference envi- 
ronment to be mitigated. In principle, a complete cross-layer solution would jointly address the routing, 
scheduling and PHY-layer issues all at once, but such a problem is very complex and remains a topic of 
future investigation. 

The rest of the paper is organized as follows. Section [TT] describes the assumed network-level model, 
introduces the definition of the network connectivity metric and formulates the general optimization 
problem to be solved. Section [III] then presents the link-level MIMO model, which includes a description 
of the time-evolution of the MIMO channels and the outdated CSI. Section [TV] describes the max- 
connectivity beamforming solution, and Section [V] puts everything together in the scheduling algorithm, 
with an appropriate power control iteration to reduce the transmit power on each link to its minimum 
possible value. Performance bounds on connectivity and sum transmit power are derived in SectionrVTl and 
results from a collection of simulation studies are presented in Section IVIII to illustrate the performance 
of the algorithm. Section IVIIII then concludes the paper and summarizes our results. The notation used 
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TABLE I 

Summary of Notations 



N 


Number of nodes in the network 


M 


Number of antennas per node 


N L 


Number of links in the network spanning tree 


Creq 


Data rate requirement for active links 


N s 


Number of time slots in a frame 


SINR k 


The signal-to-interference-plus-noise ratio of the fcth link 


SINRt 


Threshold for SINRk to guarantee the data rate requirement 


-Pout,k 


Outage probability for the fcth link 


u 


Network connectivity metric 


Xk 


Data symbol transmitted over the fcth link 


tfc 


Transmit beamformer for the fcth link 


P k 


Transmit power of the fcth link 


dik 


Distance between the transmitter of link i and the receiver of link fc 


a 


Path loss exponent 


Hit 


Channel matrix between the transmitter of link i and the receiver of link fc 


Hi fc 


Estimate of the channel between the transmitter of link i and the receiver of link fc 


Ch,l 


Largest singular value of H| fc 


Eik 


Channel perturbation for the link between the transmitter of link i and the receiver of link fc 


7 


Channel temporal correlation coefficient 


n fc 


Additive noise vector at receiver of link fc 


Wfc 


Receive beamformer at receiver of link fc 




Scheduling priority of link fc 



in the paper is summarized in Table U 

II. Network Model and Connectivity Definition 

In this section, we provide a model for the network configuration and the physical layer channel 
assumed in this paper. We also define a connectivity metric for the network, which is the performance 
objective we wish to optimize. 

A. Network Configuration and Assumptions 

We consider a multi-hop wireless network with a set of N nodes, each of which is equipped with 
M antennas (the assumption of an equal number of antennas is not strictly necessary, but simplifies the 
presentation). We assume that a source node wants to share a streaming message with all other nodes in 
the network through a pre-defined routing tree, as depicted for example in Figure Q] To avoid congestion 
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and maintain a constant average data flow from the source to all nodes, each link must achieve a certain 
minimum data rate with high probability. Performance beyond that achievable with simple TDMA-based 
protocols is possible with the availability of multiple antennas, but co-channel interference must be 
managed through appropriate scheduling, power control and transmit/receive beamforming. The ability 
of the network to achieve these goals depends heavily on the accuracy of the CSI available to the 
scheduler, as well as that available to each link of the network. 

The scheduling and transmit parameter design are centralized, and are based on outdated CSI fed 
back to the scheduler from the individual links. Consequently, we avoid the use of spatial multiplexing 
and assume that the data on each link is transmitted via a single data stream using a single transmit 
beamformer. The signal-to-interference-plus-noise ratio (SINR) of this data stream determines the rate 
of the link. Design of the receive beamformers and power control is performed by the receivers on each 
link, who are assumed to be aware of the statistics of the interference and the instantaneous effective CSI 
(channel times transmit beamformer). The source message is assumed to be broken down into packets, 
and the transmission is conducted over multiple frames. Each frame is divided into time slots, and each 
link in the routing tree is assigned one time slot in which to transmit its current packet. The scheduler 
determines the number of time slots and which links are active in each time slot, in order to optimize the 
connectivity of the network (defined below). The problem with dividing the frame into too many slots 
(emphasizing time rather than space-time multiplexing of the links) is that, for a fixed frame duration, 
the slot length is shorter, thus requiring a higher SINR to achieve the same throughput over the frame. 
This is the fundamental trade-off: fewer slots means more interference, but a lower SINR requirement; 
many slots means less interference, but a higher SINR requirement. We want to find the optimal number 
of slots for the best performance. 

B. Link Throughput Requirement 

For a given frame, we assume that each of the active links will be allocated a single time slot, and 
during this time slot, each transmitter can occupy the entire available system bandwidth to send its data 
packet to the intended receiver. Define C req to be the minimum rate at which a link can be considered 
connected, and let the number of slots in a frame be N s . For link k to meet the rate requirement, its 
SINR should satisfy the following condition: 

SINR k > SINRt, (1) 

where SINR t = 2 NsCraq — 1. Due to fading and co-channel interference, the actual SINRk may be 
smaller than SINR t . A link is said to experience an outage when the SINR at the receiver is smaller 
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than the threshold SINRt- The outage probability for link k is thus 

Pout,k = Pi{SINR k < SINR t }. (2) 

C. Definition of Network Connectivity 

We assume that the main goal of the network is robustness rather than throughput. We want to allocate 
the resources of the network (transmit power, beamformers, scheduling) such that the network remains 
connected with the highest possible probability, where the term "connected" implies that each link is able 
to communicate at or above the minimum acceptable rate C req . This in turn requires that an active link 
must have a certain minimum SINR. We define the connectivity as the likelihood that all the links in the 
network can achieve a SINR that allows transmission at or above the desired minimum data rate. Note 
that this is significantly different from simply requiring that the average rate be above some threshold, 
which is the approach taken in prior work on fast fading channels. 

If the network interference has been properly managed so that its impact on a given link is negligible, 
then the probability of a successful transmission on a given link is independent of the other links, and 
the connectivity metric can be defined as 

U = H(l-P out , k ). (3) 

fc=i 

According to this definition, the network connectivity is equal to the probability that none of the links in 
the network experiences an outage during the transmission frame. Assuming a network of Ni links and 
a frame with N s slots, SINRk can be expressed as a general function of the transmission parameters: 

SINR k = f(t h --- ,t NL ,Pi,--- ,Pn l ,N s ,S\-H), (4) 

where G c^ xl anc j {t k \k = 1, • • • , Nl} denotes the link beamformers, {Pi, • • • , Pn l } is the transmit 
power allocated to each link, Hjfc £ £_MxM an( j _ k = 1, • • • , Nl} represents the channels 

between the transmitter of link i and receiver of link k, S = {s±, • • • ,sn l } indicates the link scheduling 
scheme, s k is the slot number allocated to link k, and s& € {1, • • • , N s }. Based on (01), the connectivity 
can thus be expressed as 

Nl Nl , 

U = 1] E{I (SINR k )} = H I{SINRk\U) P {U)dU, (5) 

fc=l k=\ J 

where p(H) denotes the probability density function (PDF) of T~L and /(•) is the indicator function defined 
as 

' 1 SINRk > SINR t 
I(SINR k ) = { k ~ (6) 

SINR k < SINRt. 
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Under this definition, when the connectivity of the network is high (near one), the total network throughput 
would be approximately bounded below by the number of network links times C req . To achieve a higher 
total throughput, C req could be increased, but this would likely reduce the connectivity. A desirable 
trade-off between connectivity and throughput could conceivably be reached by appropriately adjusting 
the minimum link rate C req . 

The optimization problem associated with maximizing the connectivity of the network could thus be 
expressed as 

max U (7) 

subject to ||tfc|| = 1, 

Pk<P max , Vfce{i,--- ,N L }, 

where P max denotes an upper bound on the transmit power. Direct optimization of (O is intractable, 
since one cannot find a closed-form expression for the probability of (O in terms of the parameters of 
interest. Instead, we choose to solve the following related optimization problem: 

max V* l(SINR k \H) (8) 

subject to ||tfc|| = 1, 

P k <Pm, x , Vfc€{l,-- - ,N L }. 

The relationship between this simplified problem and the original optimization is analogous to the 
relationship between the geometric and arithmetic means. Instead of maximizing the product of the 
link connection probabilities, we maximize its sum. 

We divide the solution of ([8]) into two sub-problems: (1) scheduling and (2) beamformer design and 
power allocation. We solve sub-problem (2) for different results of sub-problem (1) to determine which 
scheduling result is best. The transmit beamforming problem is discussed in the next section, and the 
scheduling algorithm is described in Section [Vj 

III. Interference Channel Model with Outdated CSI 

A. Interference Channel Model 

Assume without loss of generality that links 1, • • • , N k are simultaneously active with link k during a 
given time slot (k > N^). The signal at link k's receiver can be expressed as 



[K ^ \Pi 

Yk = \ -i^UkktkXk + 2_^\\ ^ H ifctjXi + n fe , (9) 
V kk i=l 
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where the transmitted symbol xj. is a unit-magnitude data symbol, dik is the distance between the 
transmitter of link i and receiver of link k, a is the path loss exponent, and G C Mxl is an additive, 
spatially white noise vector with co variance given by E{ri£.n|^} = a 2 I and I is the M x M identity 
matrix. 

Assuming receiver k employs a beamformer G C Mxl , the SINR for link k is given by 

E< 



SINR k 



E< 



(10) 



where = Y2i=i J^Hjfct^Hjfctj)^ + f 2 I- Assuming the receiver knows the covariance matrix 
and the channel vector H^t^, the optimal that maximizes SINRk is given by ll27l : 



w fc = Q^Hkktk, (11) 



and the resulting SINR for link k can be expressed as 



SINR k = ^tfH&Qj^Hfcfct*. (12) 

a kk 



B. Outdated CSI 



The CSI at the scheduler will be outdated due to the time required for this information to be fed back 
from the network nodes. To quantify the CSI uncertainty due to the feedback delay, we adopt a first 
order Markov model to characterize the time variation of the channel [281 



H ifc = y/l^H-ik + ^E ik , i, k = 1, . . . , N L (13) 

where H^. denotes the channel matrix during the data transmission period, Hf fc represents the channel 
feedback, and G C is a perturbation matrix. The elements of H| fc and Ejfc are assumed to be i.i.d 
complex Gaussian random variables with distribution CM(0, 1). The coefficient 7 is used to determine 
the level of uncertainty in the CSI at the scheduler. In effect, under this model the scheduler assumes a 



Gaussian distribution for Hjfc, with mean i/l — 7H?. and independent entries with variance 7. 
Substituting (fT3l into (IT2b . the SINR at the receiver of link k can be expressed as 

SINRk = J^tf (yi^H e fcfc + ^Ekk) H Qk 1 (VT^Ktk + V7E fcfc ) t fe , (14) 

where Q fc = 4 (^ /T ^ H ^ + \/7 E ^) (v^tH^ + V^)^ + a 2 I. It is observed that 
SINRk is a function of the channel set % = {H|, \i, k = 1, . . . , iV^} and the channel perturbation set 
£ = {Ejfc|i, k = 1, . . . , A r / / }. Given the channel set SINRk is a random variable, the distribution of 
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which depends on the elements in £. The conditional expectation of SINR k with respect to £ is given by 
nSINR k \H} = l{ J£-if (v^H^ + ^E kk ) H Q k 1 (yi^H e fcfc + ^B kk ) t* j (15) 

= ^((l-7)tfH e fc fE{Q- 1 }H| fc t fe + 7 tr(E{Q- 1 })), 

where step (a) is due to the fact that perturbation matrices E ik (i, k = 1, . . . , Nl) are independent of each 
other and tr(-) denotes the trace operator. Note that the use of the expected value in (031 ) is due to the 
fact that the CSI in % may not be precisely known. According to (fT3T ). the channel may have changed 
from the time it was reported since the network nodes may be mobile. If precise CSI is available, then 
the expectation can be dropped and the instantaneous value of % can be used instead. Calculation of 
the term eIq^ 1 j is very complicated, so instead we use the following lower bound based on Jensens 
inequality |29l Lemma 4]: 

E{Q^} ^E{Q fc } _1 




where A y B denotes that A — B is a positive semidefinite matrix. In the above calculation, we use 
the fact that E^tj G C Mxl is a complex Gaussian random vector with distribution E^t,, ~ CM (0,1). 
Substituting (fl6l ) into (031 ). the conditional expectation of the SINR k is lower bounded by 

p, // Nk p \ Nk p Y 1 

E{5/Arit! fc |^}>-A(l- 7 )tfH^ 7 £^ + a 2 I + (1 - 7 ) ^ -±H| fe t,tf Hff U% k t k 

a kk \ \ i= i a ik / j =1 a ik ) 

+ 7tr ((( 7 ^ + - 2 ) I + ( 1 -7)E^H| fc t l tfH^ V (17) 
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IV. Transmit Beamforming for Connectivity 

In this section we introduce how the transmit beamformer and power are calculated for a given link 
schedule. Consider a scenario in which K links are transmitting simultaneously. In this case, the number 
of links that can meet the desired rate requirement depends on the beamformer and transmit power that 
each link adopts. Define E[{S I N R^H} as the right hand side of (fT7T) . The problem of finding the 
optimal beamformer and transmit power based on the outdated CSI % can be stated as: 

K 

^iiMSINR^H}) (18) 
subject to ||tfc|| = 1, 

Pk<Pm, x , Vfc€{l,--- ,K}. 

The indicator function in (fT8l) is not continuous and thus the problem is difficult to solve with standard 
optimization algorithms. Using the following sigmoid approximation |[30l : 

i(E l{ SINR k \U}) - 1 + e _, (M ^ } -^ t) , d9) 

where /3 is the approximation parameter, the problem can be converted to finding the maximum value of a 
constrained continuous nonlinear multivariable function. Replacing I (Ej { SINR^ \H}) with I (Ej { SINR^ \H}), 
the optimization problem in (fT8l ) can be approximated as 

K 



subject to ||tfc|| = 1, 

Pfc<iW, VfcG{l,-- - ,#}. 

Note that when /3 is small, the sigmoid function /(E^jS'/iV.Rs^l'H}) is smooth. As (3 — > oo, the sigmoid 
function approaches the indicator function. Starting with relatively small values for j3 and then increasing 
it in several steps, the solution involving the indicator function can be found. For each fixed value of /3, 
the problem in (1201 can be solved numerically using, for example, the active-set method. The algorithm 
can be initialized with arbitrary power allocations, and by setting the transmit beamformer t& equal to 
the principal singular vector of HjL. 

V. Scheduling for Maximum Connectivity 

In this section, we propose a scheduling algorithm that maximizes the connectivity metric defined 
earlier using the concept of coloring from graph theory |fl9l , ||20l . The algorithm is based on the use of 
the interference and collision graph (ICG) of the network 0, lfT8l - |[20l . 
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A. The Interference Graph and Greedy Coloring 

We use an ICG to model the relationships between the active links. In the ICG, each vertex represents 
a directional link in the transmission graph. We define two types of neighbors in the ICG. The first 
are interfering neighbors, which represent links that could be simultaneously active and hence interfere 
with one another, and the second are colliding neighbors, which represent links that cannot be active at 
the same time. In our application, colliding links include those that share the same transmitter, or those 
where a transmitter in one is a receiver in the other. All links that are not colliding are considered to 
be interfering, although the amount of interference between two given links could be low if they are far 
apart. In the paragraphs below, we more precisely define the ICG and concepts related to our scheduling 
algorithm. 

Interference and collision graph: The interference and collision graph Gi can be defined based on the 
transmission graph Gt- A given link k in Gt is represented by a vertex Vk in Gi. Suppose for link k 
that node tk is the transmitter and is the receiver, and suppose for link I that ti is the transmitter and 
H is the receiver. Links Vk and vi are colliding links if any of the following are true: tk = U,tk = ri 
or ti = rk (note that for the multicast tree we assume, = r\ will never occur). An edge between 
two vertices in Gi represents that the two corresponding links are colliding links and they could not be 
assigned to the same time slot; if two vertices in Gj do not share an edge, they represent interfering links 
which can be assigned to the same time slot, provided that the resulting interference could be managed. 
As an illustration, Fig. |2] represents a partial ICG for the network of Fig. \T\ where the colliding links and 
interfering links are connected with solid and dashed edges, respectively (for the sake of clarity, only 
edges associated with v\ and V2 are plotted; the remainder of the edges can be generated in a similar 
fashion). For example, links V2 and v\ share the same transmitter node 1, so they are colliding links in 
Fig. |2j the receiver of link i>2 is the same as the transmitter of link v&, so V2 and are also colliding 
links. 

Coloring: In our application, "coloring" refers to the process of assigning time slots to the network 
links, or equivalently, to the nodes in the interference graph. Given a set of colors in the discrete set C 
(colors can be considered as distinct non-negative integers), a coloring of the graph G is an assignment 
of the elements (or colors) in C to the vertices of G, one color for each vertex, such that no adjacent 
vertices occupy the same color. A greedy coloring enumerates the vertices in a specific order v%,...,v n 
and assigns Vk to the smallest color that is not occupied by the neighbors of Vk among ui, ... ,i>fc-i. 
The vertices can be ordered according to their edge degree, which is the number of edges incident to the 
vertex ||3T1 chap. 5]. To apply coloring to the ICG, we need to define an order for the vertices. Before 
we proceed to the scheduling order, some related definitions are necessary. 
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Scheduling freedom: For a link k, the scheduling freedom Fk is the number of available colors that 
can be allocated to this link. The higher the value of Fk, the higher the possibility that link k will be 
allocated to a "good" color (one that leads to low interference). 

Collision degree: Given a link k, the collision degree Cry,k is the number of its colliding neighbors in 
the ICG. 

Constraint and free color sets: For link k, the interfering color set Vj^ includes colors occupied by 
the neighbors that interfere with Vk- The unavailable color set T>uk includes the colors occupied by the 
neighbors that collide with Vk- The free color set is defined as T>pk = C — (Pi,k U^ , c/,fc)> and corresponds 
to the set of colors that could be assigned to link k without causing any interference (or collisions). The 
constraint color set is defined as T>c\k = T^i,k ~ (Pi,k V\^u,k)- The colors in this set can also possibly 
be assigned to link k, but with some additional interference that would have to be mitigated through 
beamforming. 

Scheduling priority: The scheduling order is determined using the largest singular values Uk,\ of the 
channel matrices H| fc . The higher the value of Cry,k or the smaller the channel gain g k = a\ Jd^. k , the 
more likely it is that link k will be affected by interference. Such a link will have fewer colors it could 
be assigned to, and hence a smaller value of Fk. To increase the likelihood that links with low scheduling 
freedom can be allocated a good color, the scheduling priority of link k is defined as 

S P ,k = C D ,k ■ W + —, (21) 

9k 

where W is a constant larger than max^ 

B. Scheduling Algorithm for Connectivity 

Based on the above definitions, we propose here a scheduling algorithm for optimizing connectivity. 
The algorithm assumes a particular value for the number of slots N s , and is repeated until a value of N s 
is found that maximizes the connectivity metric. The minimum possible value for N s is the maximum 
collision degree max^ Co,k over all vertices in the ICG. The algorithm begins by ordering the links 
according to their scheduling priority, and then assigns a color to them one-by-one, from highest to 
lowest priority. Consider the link at position m in the priority ordering. If link m can be added to a color 
that already has had other links assigned to it, and if the beamformers and power levels for these links 
can be adjusted to accommodate link m without causing any of them to drop below SINR t , then link 
m is added to this color. If there are multiple colors for which this is true, it is added to the color that 
requires the smallest increase in transmit power to accommodate it. If the addition of link m to any of 
these colors causes one of the links (including possibly link m) to drop below the SINR threshold, and 
if there exist free colors that have not had any links assigned to them, link m is assigned to one of the 
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free colors. If no free colors exist, then link m is assigned to the color that caused the smallest number 
of links to drop below the threshold (and with the smallest increase in power, in case of a tie). In this 
latter case, it is hoped that the power control algorithm described in the next section will reduce the 
interference sufficiently so that all links assigned to the color will end up being active. A more precise 
mathematical description of the algorithm is given below. 

1) Let Vfc be the right singular vector of H|, corresponding to the largest singular value, and let 
Pk = Pmas- be the initial transmit power allocated to link k. Assuming a value for N s , initialize 
the active link set A = {k \ P ma , x gk > SINR t , k = 1, . . . , Nl}. Links that do not qualify for A 
cannot meet the desired target rate for the given value of N s . Initialize the transmit beamformers 
as T = {vfc|fc E A}, the transmit powers as V = {P^ = P ma , x \k E A} and the color set 
C = {1, 2, • • • , N s }. Let Vc = {Pc,i: ■ ■ ■ i Pc,N B } denote the sum transmit power of the active 
links in each color, and let Mc = {Nc,\, • • • , Nc,n s } represent the number of links that are unable 
to meet the the target SINRt for each color. Initialize these sets to contain all zeros. The initial 
schedule S = {sk\k E A} is also set to zeros. Construct the ICG based on the relationship between 
the active links. Compute the scheduling priority Sp ; k of the vertices v k for k E A . 

2) Select the link with the highest scheduling priority: m = argmax^g^ Sp t k, and construct the free 
color set T>p iTn and the constraint color set T>c,m for link m. 

If V c , m = 4> and D F>m / 

Assign link m to color j = minj e x> Fm i, set s m = j, P m = Pqj = SINRt , and 

skip to step 5. 
else if T>c, m = 4> an d Dp, m = 4> 

There aren't enough colors to avoid collisions between the active links, so the 

algorithm must be restarted with a larger value for N s . 

else 

For i G T>c t m, construct the link set C mj i = {k \ Sk = i for k G A}, which 
contains the links currently assigned to color i. 

end. 

3) For each i G P^ m , assume the links in the set C m ^ (J m are transmitting simultaneously, and use 
the transmit beamforming algorithm of Section [IV] to find the new beamformer and transmit power 
sets T m ,i = {t l m ,k\ k e An,i U m }> Pm,i = {PJn,k\ k e ^m,i U m } ■ For lin k set C nhi U m, based on 
7~m,i, V m ,%i calculate the expected number of links that will be unable to meet the SINR threshold 
with link m added: N m ^ = Ylkec U m — l(^l{SINRk\H\)\, and calculate the updated sum 
transmit power due to the addition of link m: P m ,i = Sfce£ \jm,PL fc- ^ et ^^m,i = N m ^ — Nc,i 
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to be the number of links that will drop below the SINR threshold if link m is added to color i, 
and define AP m ,; = P m ,i — Pc,i to be the additional transmit power required to add link m to 
color i. 

4) Find j = arg minj e x) c m AN m> i. If more than one color corresponds to the minimum AN m> i, select 
the color with minimum A? mi j. 

If V F<m = 4> or AN m ji = 

Assign link m to color j' , set s m = j', N c y = N m f, Pqj' = P m ,j'- Use T m j, 
V m j' to update the components of T and V which correspond to links in the set 

An,/ U m - 
else if T>p jTn ^ 4> an d AiV m j > 

Assign link m to color j = minj g x> Fm i and set s m = j, P m = Pcj = SI ^ Rt . 

end. 

5) Set Sp t m = 0, and repeat step 2 until each vertex v^, k £ A, is allocated a color. 

Once the scheduling is complete, the active links will transmit data according to the scheduling result 
S using the beamformers in T and, at least initially, the transmit powers in V. As explained below, the 
actual transmit power for each link will be fine-tuned based on feedback from the receivers. 



C. Local Power Control for Active Links 

Since the nodes are energy limited, to extend the lifetime of the network and to reduce the mutual 
interference caused by the co-channel links, the transmit power of each link should be minimized under 
the constraint of the QoS requirement. Due to the approximation in (fT9l ), the use of the lower bound in 
(fTTT ). and the presence of outdated CSI, the actual SINRk based on V and T will not be exactly equal 
to the threshold SINR t . In most instances it will be greater than SINR t due to the use of the lower 
bound in (fTTT ). but in some rare cases it can be below the threshold. To remedy this latter situation, we 
reduce the transmit power of any links whose SINR exceeds the threshold, which reduces co-channel 
interference and the transmit power consumed by the network. For a given time slot t € C, the network 
schedule S assigns links in set £ t = {k\sk = t,k £ .4} to transmit simultaneously. The power control 
algorithm steps through each link in the time slot, reducing power for the link if its SINR exceeds the 
threshold. A given link may be revisited several times, since reductions in transmit power for other links 
reduces the overall interference, and may allow further reductions in transmit power for the link. This 
process is assumed to repeat a maximum of -/V a times. If, after all A^ a iterations, there are any links whose 
SINR is below the threshold, these links are declared to be in outage, their power is reduced to zero, 
and an additional iV& iterations are performed to reduce the transmit power even further. A mathematical 
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description of the power control algorithm for time slot t is described below. 
Iterative Power Control Algorithm 

1) Initialize the transmit power and beamformer of link k G Ct with V and T, set maximum iteration 
lengths N a and iV&. 

2) Link k's receiver (k £ C t ) calculates SINR k based on £12]). If SINR k > SINR t , the receiver 
for link k informs the transmitter to reduce P k to P k = siNRk}siNR t • 

3) 7V a = N a — 1, if N a > 0, go to step 2. 

4) If SINRk < SINR t for any A: G £ t , set = and repeat step 2 for another iVj, iterations to 
further reduce the transmit power. 

VI. Performance Bounds 

In this section, we derive two performance bounds that can be used to evaluate the limiting behavior of 
the proposed algorithms. The first is an upper bound on the network connectivity metric, and the second 
is a lower bound on the average sum transmit power. Both bounds are derived under the assumption that 
the CSI is perfect, and each active link is free of interference. When 7 is small, the connectivity bound 
should match the performance of the proposed approach if the network interference has been properly 
accounted for. This is not the case for the bound on transmit power, however, since the interference 
mitigation results in beamformers that require excess power to achieve the rate threshold. The difference 
between the required transmit power and the lower bound represents the price paid for the enhanced 
connectivity that results from operating the system as an interference network. Due to the complexity of 
calculating the bounds, expressions are derived only for the cases M = 1,2,4. 

A. Upper Bound on Network Connectivity 

The connectivity bound is derived assuming the absence of interference, and perfect CSI at the scheduler 
(7 = 0). Each transmitter uses maximum power and selects the principal singular vector of the channel 
matrix as its beamformer. When the receiver is free of co-channel interference, the resulting signal-to- 
noise-ratio (SNR) for link k is given by 

P n 2 

If we define P t k = PT{SNRk < SINRt}, then the upper bound on connectivity can be expressed as 

U^ 1 = U(l-Pout,k)>U. (23) 
fc=i 
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The squared singular value a\ x corresponds to the largest eigenvalue of the central Wishart matrix 
H fcfc Hf fc . Define X k . i — <Jf.i, and note that the cumulative density function (CDF) of A^i is given by 
(321 eq. (6)]: 

M^^ /rv (24) 

(n",r 2 (i)) 

where we have dropped the subscript k since the distribution is assumed to be identical for each link, 
T(-) is the gamma function, 3?(A) is an M x M matrix defined by 3?(A)j j = j(i + j — 1, A), and the 
lower incomplete gamma function 7(n, x) has the following series expansion: 



V fc=o / 



7(n,A) = (n- !.)!( I. - > ^ —<~ A ) . (25) 

For M = 2, ([24]) reduces to: 

Pr{A a < A} - det W A ) 2 * 2 



n; =1 r 2 (j) 

7(1,A) 7 (3,A)- 7 2 (2,A) 

r 2 (i)r 2 (2) 

= l-e" A (A 2 + 2) + e" 2A . (26) 
Substituting (l26l ) into (l23T ). the connectivity upper bound can then be calculated as: 

U%= 2 = e- ^ A — II ~ e~ X - k + 2 ) > ( 27 ) 

fc=l 

where A m i nj fc = SIN p td >*>* a ' i s the minimum value of the channel gain that can guarantee SNR^ > 
SINRt. 

For M = 4, after some cumbersome algebra, the CDF of the largest eigenvalue Ai can be expressed 
as: 

det($(A) 4X 4) 



Pr{Ai < A} 



nj =1 r 2 (i) 

l-/i(A)-/ 2 (A)-/ 3 (A)-/ 4 (A), (28) 



where /i(A), /2(A), /3(A), /4(A) are defined in the appendix. The network connectivity upper bound is 
then computed as 

N L 

Ug~ A = Y\ (/l(Aniin,fc) + /2(Amin,A:) + /3(A m in,fc) + /4(Amin,fc)) • (29) 
k=l 

For the single antenna case (M = 1), the interference-free SNR at the receiver is given by 

SNR k = Pm ^ h f l \ (30) 
a kk a 
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where \hkk\ is a Rayleigh random variable. Define h = \hkk\, so that Pr{/i < A} = 1 — e A2 . The 
probability of a successful transmission for the link can be expressed as 

Pr{SNR k > SINRt} = 1 - Pr jfc < a/A~^| = e -W, (31) 

and thus the network connectivity is given by 

^=i = e -£^W. (32) 



Comparing (1271 ) and (1291 ) with (1321 . it can be observed that 

jjM=2 Nl 

R2 = = J] A min, fc " e~ W +2 , (33) 

Ra = T f M= i = Y\ eAmm ' fc (/l(A m in.A:) + /2(A m i n> fc) + /3(A m in,/c) + /4(A mirl .fc) J . (34) 

It is easy to verify that both R2 and R4 are monotonically increasing functions of X m m : k and larger than 
1. i?2 and i?4 represent the connectivity gain provided by the use of multiple antennas. Since X m - m ^ is 
proportional to d^ k , the larger the link distance, the greater the connectivity gain offered by the use of 
multiple antennas. 



B. Lower Bound on Average Sum Transmit Power 

A lower bound for the average sum transmit power of the network can be obtained under the assumption 
that each active link selects its beamformer as the right singular vector of the channel with largest singular 
value, and assuming there is no co-channel interference between the links. 

For link k, given Ai, the power allocation at the transmitter is: 

pM = \ A, > A ^W (35) 

I 0) Ai < A m i ni fc 

The average transmit power of link k can be obtained by averaging P^ 1 over the random variable Ai. 
Denote the probability density function (PDF) of Ai as f\ 1 (\), then f\ 1 (X) can be explicitly obtained 
by taking the derivative of (l26l ) or (l28l) with respect to A. 
When M = 2, f^ =2 (X) is given by 

f^= 2 (X) = e" A (A 2 + 2) - 2Ae" A - 2e^ 2A . (36) 
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The average transmit power for link k is calculated a 

M 2, _ /'' V: ~ SINR t d kk a2 f M 2, 



E{Pf = /. I fcfc /jT 2 (A)dA 



SINR t dl k a 2 [ Ae~ A - 2e" A + ^ (e~ A - e~ 2A ) dA . (37) 
Jx mi „. k A 



rA„ 
'A mln , fc 

In (I37T ). the integration of the first two terms can easily be found. To calculate the term J ^y-dA for 
constant c, we expand the exponential function /(A) = e _cA using a Taylor series. The Taylor series at 
the point A is 

JE ( — r\ n P~ cX 

*"(*) = E - n , (A - A) n . (38) 

n=0 

With the help of ([38]), the term / ^dA in OH) can be further expressed as 



n=0 fc=0 

r + E iz£ ^ f ^ + E it) (-A)'-^*-') ,» 

n=l ' * \ k=l v 7 / 

°° f „\n a — cX ( n 

cA 



lnA +E (<-*>" b- a + E (I) <-*>"-" T • <3« 

n=l \ fc=l V 7 > 

and the average transmit power for link k can be written using the following expression.:: 

E{Pf = 2 } = SINR t d% k a 2 {{\ - A max )e- A — - (1 - A min , fc )e- A — fc + 2 5l (l, A max ) 

-201 (1, \ m in,k) ~ 01 (2, A max ) +01 (2, \ m in,k))- (40) 

Similarly, when M = 4, the PDF of Ai can be expressed as: 

/a M=4 (A) = /i(A) + / 2 (A) + / 3 (A) + /1(A), (41) 

where the definitions of /{(A), /{(A), /i(A), /{(A) can be found in the appendix. The average transmit 
power for link k in this case is given by 

,m=4t _ /" Amax SINR t d% k a 2 M=A 



= SINR t d a kk a 2 ^ maX I(/;(A) + / 2 (A) + / 3 (A) + /;(A))dA. (42) 

'in principle, A max — > oo in this equation, but to numerically evaluate the integral, we simply choose a large enough value 
such that the integrand is essentially zero. 

2 For gi(c, A), faster convergence of the series can be obtained if A is selected as A = A ™ ax . 
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Based on (I47l)-(l50t in the appendix, the average transmit power in (1421 ) can be evaluated as 

^{Pk I=4 } = SINRtd^a^hiX^) - /i(A min , fc ) + / 2 (A max ) - / 2 (A min , fc ) + / 3 (A max ) 

— /3(A m in,fc) + /4(A m ax) — /4(A m i nj fe) ) . (43) 



For the single-antenna case, we average the transmit power over the channel gain h^f. to determine 
the average transmit power of link k as 

K{P k M =1 } = T maX 2SINR f kka2 e~*dh 



= 2SINR t d a kk a 2 (s~(A max ) - «), (44) 
where the function s(x) is defined as 



six 



„-x 2 00 „2n 



cfc = In* + £(-!)»_. (45) 



x / -— ' 2n ■ nl 



Thus, the lower bound on the average sum transmit power of the Nl links is given by 



s k=i 



VII. Simulation Results 



(46) 



For our simulations, we consider a network with N = 30 nodes, uniformly distributed in a 25m x 25m 
area as shown in Fig. [3] The node represented by a square is the source node, and the edges represent the 
Nl = 29 links in the multicast tree. Also in Fig. [3j the scheduling result for a single channel realization is 
provided when N s = 3, M = 4, 7 = 0.04 and C req = 0.9bps/Hz; the links that have the same color have 
been scheduled to transmit in the same time slot. We assume noise with unit power, a maximum transmit 
power of -P max = 25dB for each node, and a path loss exponent of a = 2. The connectivity performance 
and the sum transmit power of the network are averaged over 300 independent channel realizations and 
the performance for different M, N s and 7 is provided. Note that although the highest collision degree 
for any of the nodes is 4, the minimum number of slots considered for a given frame is 3, which means 
our scheduling algorithm can use only 3 time slots to completely avoid the link collisions. 

In the power control algorithm, only the first few iterations play an important role in the algorithm 
performance. This is illustrated in the example of Fig. [4] for a case with M = 4, C req = 0.9 bps/Hz, 
7 = 0.1 and N s = 5. We see that all of the power reduction occurs for {A^, N^} < 3. In the simulation 
results that follow, we set N a = 3 and = 2. The minimum rate requirement C req in the following 
simulations is assumed to be adjusted to take into account the overhead due to channel estimation and 
feedback to the source. The QoS requirements are set to be C req = 0.9bps/Hz for M = 4, C req = 
0.5bps/Hz for M = 2 and C req = O.lbps/Hz for M = 1, and plots for both 7 = 0.01 and 7 = 0.04 are 
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included. We show results for different C req with each M, due to the strong impact of the number of 
antennas on performance. When the connectivity probability is near 1, the approximate total throughput 
of the network can be found by multiplying C req by Nl — 1 and the connectivity probability. 

We compare the performance of three cases: outdated CSI (OCSI), local perfect CSI (LPCSI) and 
global perfect CSI (GPCSI). For OCSI, the source node uses the outdated global CSI % as the input to 
the scheduling algorithm, and the links transmit according to the scheduling results S, T and V provided 
by the source. In LPCSI, the CSI of the links transmitting in the same time slot are assumed to be known 
by the other active links and, based on this local CSI, the transmit beamforming algorithm is used to 
re-optimize the beamformers for that time slot. The performance gain of LPCSI over OCSI represents 
the advantage provided by the use of local instantaneous CSI. For GPCSI, we assume the source node 
has perfect knowledge of T~L, which amounts to assuming 7 = 0. For this case, the lower bound on 
SINR in (fT8T ) is replaced with the exact SINR expression in (fl2l ). The results for GPCSI indicate the best 
performance that the scheduling algorithm can achieve, and we see that as predicted, its performance 
matches the bound in Section IVI-A1 provided that N s is chosen large enough so that the interference can 
be properly mitigated. 

Fig. [5] provides the connectivity performance and the average sum transmit power for M = 4. The 
behavior of the connectivity metric can be explained as follows. As the number of colors (slots) increases, 
the number of links that transmit simultaneously will be reduced, the interference between the links will 
thus be reduced, and the SINR at the receiver of each link will be improved. However, to guarantee the 
spectral efficiency C req , the threshold SINR t must also increase due to the shorter time slot. When the 
benefit brought by the increase in the number of slots is larger than the penalty caused by the increased 
SINR t , the connectivity of the network will increase; otherwise, the connectivity will decrease. This 
trade-off results in an optimal value for N s for each case considered. 

The behavior of the transmit power curves is slightly different, since links that cannot meet the desired 
SINR are not allowed to transmit. This obviously reduces the connectivity metric, but it also reduces the 
total transmit power. That explains why, for example, the transmit power required by GPCSI is always 
higher than that for the other cases; since it achieves a higher connectivity, more links are active and more 
transmit power is consumed. Note also that the analytical values for both the connectivity and transmit 
power bounds match those obtained in the simulation. It can be observed that when 7 = 0.04, full 
connectivity is achieved at C req = 0.9bps/Hz with a required transmit power of about 150. Comparing 
the performance of GPCSI and OCSI for 7 = 0.04, we see no impact of the outdated CSI on the network 
connectivity for N s < 5. 

In Fig. we see that the two-antenna network is able to achieve a connectivity of about 0.7 at 
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Creq = 0.5bps/Hz, requiring a transmit power of approximately 300. Compared with GPCSI for 7 = 0.04 
and N s = 5, the network connectivity of OCSI is reduced by about 25%. 

Fig. [7] presents the performance results for M = 1 with C req = O.lbps/Hz. Peak connectivity occurs at 
N s = 5, with a required transmit power of about 140, and a connectivity of only about 0.15 is achieved. 
When 7 = 0.04, we see that the connectivity performance of OCSI is reduced by a factor of three with 
respect to GPCSI. 

Comparing the connectivity performance in Figs. [5]|7J it can be observed that the optimal value for 
N s decreases with M, even though we have increased the desired throughput with M as well. The 
benefit of having nodes equipped with multiple antennas is clearly evident in terms of connectivity and 
total network throughput. Note that the benefit of nodes with multiple antennas also manifests itself in 
terms of the transmit power required to achieve maximum connectivity. Instead of a throughput gain 
of four with the same level of reliability, which might be expected when comparing the M = 4 and 
M = 1 networks, we see that there is a multiplicative benefit that results from using multiple antennas for 
improved connectivity. The additional spatial degrees of freedom reduce interference between links within 
the same time slot, they reduce the total number of required time slots, and they reduce the amount of 
transmit power (and hence interference to other links) to achieve a given throughput. The combined gain 
of these effects can only be observed by considering the joint optimization of the transmit beamformers, 
scheduling and power control. 

Further evidence of the multiplicative benefit of multiple antennas is the fact that we observe a decrease 
in sensitivity to imprecise CSI as the number of antennas at each node increases. Note also that using 
locally accurate CSI to adjust the transmit beamformers (LPCSI) has only a slight benefit in improving 
performance relative to using the beamformers based on outdated CSI (OCSI). This is due to the fact 
that only a single data stream is transmitted on each link. If multiple data streams were allowed, the 
importance of locally accurate CSI would be more critical; one could potentially improve connectivity, 
but the performance gain would be much more sensitive to imprecise CSI. 

In Fig. [8j we show the connectivity results for all three antenna sizes for the same value of C reg = 
0.6bps/Hz. The M = 4 network achieves full connectivity in this case, while the M = 1 network is 
nearly disconnected. The two-antenna case falls somewhere in between, with performance depending on 
how accurate the CSI is at the scheduler. 

In Fig. |9j the transmit power results are provided. The difference in transmit power is especially evident 
in this example. The M = 4 network requires a factor of seven times less power to achieve a connectivity 
that is five times higher than the M = 2 case. While the single-antenna network is almost disconnected, 
the fact that the average sum transmit power is non-zero indicates that at least some of the links are 
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active. 

Finally, in Figs. [T0lfT2l we compare the performance of our connectivity-based approach with one 
that attempts to choose the network parameters in order to maximize the throughput of the network. 
The optimization algorithm in this case is similar to the scheduling algorithm for connectivity, with the 
exception that the objective function in ([8]) is replaced with X)fciil°§2(l + SINR^). In Fig. [TOj for 
M = 2, we see that attempting to maximize the throughput results in zero connectivity for the network, 
while using our approach, we achieve a connectivity in excess of 0.95 for the optimal value of N s = 5. 

In Fig. [TT] the overall network throughput of the max-throughput approach is about 50% higher than 
the connectivity optimization method. This indicates that a subset of the nodes is able to communicate 
at a higher rate, but the network is disconnected and the message from the source is not reaching all of 
the network nodes. 

In Fig. [I2j the results show that the proposed approach achieves the maximum connectivity with about 
five times less transmit power than the max-throughput approach. The connectivity of the max-throughput 
approach improves when M = 4, but with an increase in required power consumption of nearly a factor 
of eight. 

VIII. Conclusions 

In this paper we investigated the use and benefit of multiple antennas in an ad hoc network with a source 
streaming data to all nodes via a multi-hop tree. Based on a novel definition of network connectivity, 
a scheduling algorithm was developed that takes advantage of the interference mitigation capabilities 
of the MIMO nodes. A key component of the algorithm is the design of transmit beamformers for 
simultaneously active nodes that optimizes the connectivity metric. Ultimately, the scheduling algorithm 
acts to break down the full network into a set of smaller interference networks whose links are able 
to be simultaneously active due to the interference mitigation provided by the multiple antennas. We 
also derived performance bounds on the network connectivity and average required sum transmit power, 
assuming zero interference and perfect CSI. These bounds represent ultimate limits for the performance 
of the network, and were used in the simulations to compare against the actual behavior of the network. 
Our simulation results indicate the significant advantage provided by multiple antennas in the ad hoc 
network, in terms of connectivity, throughput, reduced transmit power and resilience to outdated CSI. 
By considering the joint problem of transmit beamformer design, scheduling and power control, we 
observe a multiplicative benefit to the use of multiple antennas for improving network reliability. Such 
an observation could not be made by solving each of these problems in isolation from the others. 
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Appendix 



The functions /i(A), /2(A), /3(A), /4(A) in (|28T ) are denned as follows: 



/i(A) 


= (e" A + e" 3A ) 


/2(A) 


= (e -A _ e -3A } 


/ 3 (A) 


= -e' 2X (l2\ 2 


/4(A) 


= 4e" A + 4e" 3A 



6 36 



-4A 3 - ^A 5 
3 



3 A + 9 A +144' 

6e -2A_ e -4A^ 



Based on /i(A), /2(A), /3(A), /4(A), the functions f[(\), / 2 ( A), / 3 ( A), /4(A) in (HU are defined as: 
/i(A) = = (e- A + e- 3A ) (-12A - |a 3 - I A 5 ) - (e" A + 3e~ 3A ) (-6A 2 - ^A 4 - ^A c 

#2(A) _ , A q„-3A\/,i\3 , 1\5\. , /„-A „-3A\ / io\2 , 5 \4 



/ 2 (A) = = (-e" A + 3e- 3A ) (^4A 3 + -A 5 J + (e" A - e~ 3A ) (^12A 2 + -A 

/3 ( A ) = = e -2A ( 24A + 8 A 3 + f_ 5 1 a tA _ 2e -2A / 12A 2 2 4 2 6 1 , 

J3K J dX \ 3 3 18 J V 3 9 144 

/ 4 (A) = = 4e~ A + 12e" 3A - 12e" 2A - 4e" 4A . 

Define ^(c, = f x n e cx dx = J27=o c'+^n-iy. xU ~ %eCX ^ then it can be shown that 

AM 

= - 12( 52 (-1, 0, A)+0j(-3, 0, A)) -^( 52 (-l, 2, A)+ 92 (-3, 2, A) 
-i( 52 (-l,4,A)+ 52 (-3,4,A)) +6(«? 2 (-l,l,A)+3 52 (-3,l,A) 



+ y ( 52 (-l, 3, A) + 3g 2 (-3, 3, A)) + ( 52 (-l, 5, A) + 3<? 2 (-3, 5, A)) , (47) 



/ 2 (A) = / ^dA 



A 

4(3 52 (-3, 2, A) - 52 (-l, 2, A)) +i (3<7 2 (-3, 4, A) -«&(-!, 4, A) 



+ 12( 52 (-l,l,A)- 92 (-3,l,A)) +^( 52 (-l,3,A)- 52 (-3,3,A)), 



(48) 



/3(A) = / &^dX 



=24 52 (-2, 0, A) + |? 2 (-2, 2, A) + -#,(-2, 4, A) + -£fla(-2, 6, A) 

- 24 52 (-2, 1, A) - |s2(-2, 3, A) - i? 2 (-2, 5, A) - ^<fc(-2, 7, A), (49) 

/ 4 (A)= /" ^dA = 4( 5l (l,A)-gi(4,A)) + 12(^(3, A) - 51 (2, A)). (50) 
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Fig. 5. Connectivity performance and average sum transmit power of the network for different number of time slots per frame, 
with M = 4, C req = 0.9bps/Hz. 
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Fig. 6. Connectivity performance and average sum transmit power of the network for different number of time slots per frame, 
with M = 2, C req = 0.5bps/Hz. 
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Fig. 7. Connectivity performance and average sum transmit power of the network for different number of time slots per frame, 
with M = 1, Cre q = O.lbps/Hz. 
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Fig. 8. Connectivity comparison between different antenna configurations M = 4, M = 2, M = 1, with C re q = 0.6bps/Hz, 
7 = 0.01 or 7 = 0.04. 
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Fig. 9. Average sum transmit power comparison between different antenna configurations M = 4, M = 2, M = 1, with 
C req = 0.6bps/Hz, 7 = 0.01 or 7 = 0.04. 
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Fig. 10. Connectivity performance of the network for different scheduling methods, with Creq = 0.5bps/Hz for M = 2 and 
C reg = 0.9bps/Hz for M = 4. 
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Fig. 11. Average sum throughput comparison between different scheduling methods, with C req = 0.5bps/Hz for M = 2 and 



0.9bps/Hz for M = 4. 



3500 r 



300 



-W-GPCSI, M=2 
-©- Max Throughput, M=2 
-B-GPCSI, M=4 

Max Trhoughput, M=4 



CD 

O 
CL 

E 
w 



2500 



200^ 



w 1500 

CD 
CO 
CD 

< 1000 



500 



0. 



o 



0~ 



a - 
V 



-v-- 



" ~ -G- - - 



--o 

- V 



5 6 7 

Number of Slots, N 



Fig. 12. Average sum transmit power comparison between different scheduling methods, with Creq ~ 0.5bps/Hz for M — 2 
and Creq = 0.9bps/Hz for A/ = 4. 
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